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We investigate the Josephson coupling between two singlet superconductors separated by a half- 
metallic magnet. The mechanism behind the coupling is provided by the rotation of the quasiparticle 
spin in the superconductor during reflection events at the interface with the half metal. Spin ro- 
tation induces triplet correlations in the superconductor which, in the presence of surface spin-flip 
scattering, result in an indirect Josephson effect between the superconductors. We present a theory 
appropriate for studying this phenomenon and calculate physical properties for a superconduc- 
tor/half metal/superconductor (S/HM/S) heterostructure. 
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Introduction: The interplay between superconductivity 
and spin-polarized materials has potential applications in 
the emerging field of spin electronics. For this purpose, a 
high degree of spin polarization of the materials in con- 
tact with superconducting regions is desirable. The re- 
cently discovered half metals are ideal materials in this 
respect [l[ . In half metals electronic bands exhibit insu- 
lating behavior for one spin direction and metallic behav- 
ior for the other. They are thus completely spin-polarized 
magnets. Half-metallic behavior has been found exper- 
imentally in the manganese perovskite Lao.7Sro.3Mn03 
0, El and in Cr02 0] . The perovskite is particularly in- 
teresting because of its ability to form high-quality het- 
erostructures with high-T c cuprate superconductors 0. 

The superconducting proximity effect in spin-polarized 
materials has attracted considerable attention recently in 
the context of superconductor/ ferromagnet heterostruc- 
tures 0, S H IS [H3 d ■ The singlet pairing amplitude 
shows oscillations with a wave vector matching the spin 
splitting of the Fermi wave vectors in the ferromagnet 
|lllll2| . The magnitude of this proximity effect decreases 
with increasing spin polarization. In the extreme case of 
a completely spin-polarized material the singlet proxim- 
ity effect is absent. Consequently, the Josephson cur- 
rent between two singlet superconductors separated by a 
half metal is expected to be exactly zero. In this Letter 
we show that this is not necessarily the case. We pro- 
pose a mechanism which leads to a nonvanishing S/HM/S 
Josephson effect. 

The indirect Josephson effect requires the interplay of 
two separate interface effects: spin mixing (or spin ro- 
tation) and spin-flip scattering. The former, represented 
by the phase difference 8 between waves of opposite spin 
orientations reflected from a spin-active interface, intro- 
duces triplet correlations at the superconducting side of 
the S/HM boundary. The latter mediates these corre- 
lations to the half-metallic side. To illustrate the spin- 
mixing effect, consider the reflection of two quasiparti- 
cles, \\)k and \l)k, from a half-metallic material (which 
defines the spin quantization axis). The reflected ampli- 



tudes for opposite s pins differ in phase, |T)-fc = e* 6 */ 2 |T)fc 7 
\l)-k = e~ l6, / 2 ||)fc .13]. In a superconductor, incoming 
quasiparticles (fc) near the interface form pairs with out- 
going quasiparticles (-Ar). As |T)fe|j)-k - |J.)fe|t)-fc trans- 
forms under reflection into e te \1)k\l)-k — e _ie |j)fc|T)-fc, 
pairing states near such interfaces are singlet-triplet mix- 
tures. This property of spin mixing is intrinsic to any 
spin-active interface. If, additionally, spin-flip scattering 
is present at the S/HM interface, the resulting triplet am- 
plitudes induce equal-spin pairing correlations in the half 
metal, leading to an S/HM/S Josephson effect. Spin-flip 
scattering is expected to be enhanced e.g. due to local 
variations of the spin quantization axis or diffusion 
of magnetic moments. The importance of these processes 
was pointed out by recent experiments flij ] . 

The indirect proximity effect introduced above can also 
be relevant for strong ferromagnets. In the conventional 
description, the dispersions for spin-up and spin-down 
bands in ferromagnets are assumed to be identical apart 
from an energy splitting, given by an effective exchange 
field h 0, 0- The range of the spin-singlet proximity 
effect is drastically reduced by a strong exchange field. 
In contrast, no such suppression occurs in the case of the 
indirect proximity effect. 

Theory: Our treatment is based on the quasiclassical 
theory of superconductivity |l5l |. This theory is formu- 
lated in terms of Green's functions (propagators) which 
are matrices in Nambu-Gor'kov particle-hole space and 
in spin space. The quasiclassical propagator, g(k,R, e) 
depends on energy e, position R, and the direction k of 
the momentum on the Fermi surface. Its particle-hole 
diagonal and off-diagonal elements are denoted by spin 
matrices g and /. The quasiclassical propagator satisfies 
the Eilenberger equation [ljj| 
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with the Fermi velocity, Vf (k), and the singlet order pa- 
rameter A(R). It is essential for our purpose to deter- 
mine the spatial variation of the order parameter near 
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the interface region in accordance with the triplet corre- 
lations, which decay into the superconducting region on 
the coherence length scale. In order to ensure current 
conservation in the whole system we obtain the spatial 
variation of A(R) self-consistently, 



^ (k ou t kin ) — S x S 



A(R) = A 



(2) 



The coupling constant A and the cut-off energy e c are 
eliminated in favor of the transition temperature T c in 
the usual manner. The quasiclassical Green's functions 
are normalized according to g 2 = — ir 2 l [l5j . 

Boundary conditions: A standard method to treat 
boundary conditions for spin active interfaces is a scat- 
tering matrix formulation fill llij. However, for the 
present problem, where the number of spin channels on 
one side of the interface differs from that on the other, 
it would be necessary to use the formulation by Millis et 
al. which is rather involved. For this reason we p ro- 
ceed with an alternative but equivalent approach j^j. It 
allows us to derive explicit quasiclassical boundary condi- 
tions in terms of an auxiliary Green's function, g° , which 
solves the boundary condition for an impenetrable inter- 
face and is easy to obtain. The impenetrable interface is 
characterized by two surface scattering matrices, S and 
5, on either side of the interface. The resulting prop- 
agators on the two sides are denoted by g° and g°, re- 
spectively. At the boundary, incoming propagators, gf n , 
are connected with outgoing ones, t , via the surface 
scattering matrices by g° ut = Sg" n S^ 13]. Particle con- 
servation requires unitarity, = S~ 1 . We include the 
transmission processes through the interface via an effec- 
tive hopping amplitude f in a t-matrix approximation. 
We assume translational invariance in the plane of the 
interface. The quasiclassical hopping amplitudes from 
left to right differ in general for incoming and outgoing 
quasiparticles. However, the requirement of current con- 
servation leads to relations between these elements as 
shown in Fig. ^ 

The quasiclassical t-matrix equations read 



tin — T S-out^ ( 9in^in) j t Q ut — Sti n S , 
Lut = ^g° in T[l lj°L) , 4 
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On each side of the interface, the t matrix describes the 
modifications of the quasiclassical propagators due to vir- 
tual hopping processes to the opposite side. Finally, we 
express the full propagator in terms of the decoupled so- 
lution g°, leading to the boundary conditions for incom- 
ing and outgoing propagators, 

(j™ = g"„ + [gi t + ini} i m {<?°„ - wrij, (4a) 

{font - ™i] Lut {9out + mtI}. (4b) 
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FIG. 1: Scattering geometry illustrating the scattering chan- 
nels and the corresponding transfer amplitudes for the model 
discussed in the text. 



and similarly for g in and g |211l22l ]. In the appropriate 
limiting cases these boundary conditions reduce to those 
published previously [3 [H H H HJ . 

For reference, we also present the corresponding full 
scattering matrix which would enter the boundary con- 
ditions of Ref. Without loss of generality it can be 
written in the form 



S21 S22 
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with the transmission matrix d = (1 + 7r 2 f f t) 1 2ttt 
the reflection matrices on cither side of the interface, f 



and 



(1 +7T 2 fft )~ 1 (I- 7T 2 ffT), f = (I 



(1 
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The particle-hole structures of the surface scattering 
matrix and the hopping amplitude are S =diag[S*, S] and 

f =diag[r, §W*S_ ]. The above equations are for general 
spin structures. In the following, r is a 2x1 spin matrix, 
S a 2x2 spin matrix, and S_ a spin scalar. 

S/HM/S structure: We study a heterostructure con- 
sisting of a half metal, — Lh < x < Lh, between two su- 
perconductors, —L<x< — Lh and Lh <x<L. We inves- 
tigate the equilibrium supercurrent due to a phase differ- 
ence 4> between the superconductors, A(L) = A(—L)e 1 ^. 

As mentioned above, band splitting in the interface 
region results in a relative spin phase for quasiparticles 
with spin along the quantization axis of the half metal 
(for quasiparticles with spin in the perpendicular plane 
the corresponding effect is a spin rotation around the 
quantization axis) |l3l | . This effect can be described by a 
scattering matrix S — exp(i9a z / 2)1 at the superconduct- 
ing side of the interface, where 9 defines a spin-rotation 
angle and a z denotes the Pauli spin matrix 0,0. Gen- 
erally, the value of 9 depends on the angle of impact, ip 
[l3j | and can approach values of the order of ir for strong 
band splitting |23j. For definiteness, we present results 
for 9 = 0.757T cosip. On the half-metallic side the scat- 
tering matrix has no spin structure, £7 = 1. 

The i-matrix equations are parameterized by the hop- 
ping amplitude f and the scattering matrices S, S_, which 
are the phenomenological parameters characterizing the 
interface in our theory. We use r = (1 + S^)t cos tp, 
where tq = (r^ ,T||) T is determined by the two spin 
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FIG. 2: Self-consistent order parameter and triplet correla- 
tions in an S/HM/S heterostructure for a zero junction and 
a 7r junction. The relative signs of the pairing correlations in 
the s-wave singlet and three p-wave triplet channels are indi- 
cated. A zero junction for the singlet order parameter leads to 
a relative phase difference of n for the triplet correlations, and 
vice versa. The calculations are for temperature T = 0.05T C , 
and for t^/t-\-\ = 0.7. 
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FIG. 3: Critical Josephson current density as a function of 
temperature for an S/HM/S heterostructure. The two curves 
are for Tlj/tjj — 0.1 (dashed) and t^/t^ = 0.7 (full lines). 
The inset shows the current-phase relationships for rj.f /Tfj = 
0.7 for temperatures T/T c = 0.05 (dashed), 0.2,0.3,0.4,0.5 
(full lines from bottom to top). The unit is the Landau critical 
current density Jl = evjNfAo, with the zero temperature 
bulk superconducting gap Ao = 1.76T C . 



scattering channels from the superconductor to the half- 
metallic spin-up band. With this choice the spin rotation 
during transmission is half of the spin rotation during re- 
flection. The cos ip factor accounts for the reduced trans- 
mission at large impact angles. We present results for 
T iT /r TT = 0.7 and 0.1, 2ttt tt = 1.0, 2L H = 3£ (with 
the coherence length £ = Vf/^Tc), L 3> Lh, and cylin- 
drical Fermi surfaces (calculations using spherical Fermi 
surfaces lead to similar results). We iterate Eqs. (1) and 
(2) until self-consistency is achieved, with the boundary 
conditions (3) and (4) at the two interfaces. All our cal- 
culations are in the clean limit. 

In Fig. |2 we present the spatial modulation of the sin- 
glet order parameter and the triplet pairing correlations 
for an S/HM/S heterostructure. We compare results for 
a zero junction (t/> — 0) and a tt junction (c/> = tt). The 
spin-rotation effect at the superconducting side of the 
interface leads to local triplet correlations in the super- 
conductor of the form + f^. We quantify the triplet 
pairing correlations by the integral 

^- (!?(£)/(£, s,e)> t tanh (^), (6) 

where ?y(k) projects out the p-wave pairing amplitude, 
and is equal to the cosine of the angle between k and 
the surface normal. Spin-flip scattering induces a 
amplitude in the half metal, and leads to both and 
Fn amplitudes in the superconductor. The correlations 
are shown in Fig. for all three spin-triplet channels. 
Triplet correlations extend into the superconductor up 
to a few coherence lengths from the interface, leading to 



a suppression of the singlet order parameter near the in- 
terface. We also show schematically the s and p orbitals 
for a zero junction and a tt junction. The alignment of 
the p orbitals is determined by the direction of the sur- 
face normal. As a consequence, the relative sign between 
the p orbitals is opposite to that of the s orbitals. As will 
be shown below, this leads to a reversal of the current di- 
rection from that expected for a superconductor/normal 
metal/superconductor junction. 

We now turn to the half-metallic region in Fig. [5J The 
spatial distribution of the proximity- induced ff j ampli- 
tude shows a sign change at x = in the case of zero 
junction, but not for a tt junction. As a result, the tt 
junction is expected to be more stable than the zero junc- 
tion. Indeed, our numerical calculations show that the tt 
junction corresponds to the free-energy minimum for all 
temperatures. The equal-spin correlations decay slowly 
into the half metal, e.g. F^ (x = 0) cx 1/Lh in the tt 
junction. This behavior is similar to that observed in 
normal metal/superconductor structures. 

In Fig. 0we show the Josephson critical current as a 
function of temperature. The current density, 

J= f°° de(ev f (k) Ni(k,e))tn f (e), (7) 

J — oo 

is expressed in terms of the angle-resolved density of 
states at the Fermi surface in the half metal, iVj 

Y ■ Im(<7i"f)/7r, the electronic charge e, and the Fermi 
distribution function rif . In the inset of Fig. |21 we show 
the current-phase relationship for different temperatures. 
The current is negative for a positive phase difference <fi. 
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FIG. 4: Density of states at T = 0.05T C for quasiparticles 
with normal impact at the half-metallic side of the left inter- 
face (x = —Lh), for a) spin flip rate Tit/Tit = 0.1 and phase 
difference (j> — 0.5ir, and b) t^/t^ = 0.7 and <j> = 0.2n. The 
corresponding Josephson currents are close to the critical val- 
ues. Shown are both states carrying current in positive (full 
lines) and negative directions (dashed lines). 



For each temperature we determine the critical current 
from the maximum current magnitude in the current- 
phase relationship. The critical current has a (1 — T/T c ) 2 
dependence near T c . This is a consequence of the fact 
that the order parameter at the interface varies linearly 
with 1 - T/T c , in contrast to the bulk (1 - T/T,,) 1 / 2 
behavior. At low temperatures the critical current passes 
through a maximum and then decreases again. This 
anomaly is due to the interplay between current-carrying 
states, as we proceed to explain. 

We discuss the different contributions to the Joseph- 
son current coming from the spectral features in the 
momentum-resolved density of states N-\ in the half 
metal. The total current through the interface is dom- 
inated by quasiparticle trajectories parallel to the sur- 
face normal. In Fig. 0]we compare the spectrum of such 
quasiparticles for incoming and outgoing momenta at the 
half-metallic side of the left interface. We present results 
for T||/r|| = 0.1 and t^/t^ = 0.7. In both cases there 
is a continuum around the chemical potential (e = 0). 
On either side of this continuum there is a gap, followed 
by either additional continuum branches, or by Andreev 
bound states. The Andreev bound states in Fig. are 
closely related to the surface Andreev states discussed 
in Refs. 0, HjJ According to Eq. Q), the current is 
obtained by multiplying the curves in Fig. 0] with the 
Fermi function. At not too low temperatures the Joseph- 
son current is dominated by the negative-energy features 
below the continuum at the chemical potential. These 
features carry current in negative direction, explaining 
the negative sign of the Josephson current for positive 
phase difference. Below a certain temperature, the cor- 
responding states are fully populated, and the tempera- 
ture dependence of the Josephson current is dominated 
by the low-energy continuum around the chemical poten- 
tial. The current carried by this low-energy band is pos- 
itive and increases with decreasing temperature, leading 



to the decrease of the magnitude of the critical Josephson 
current at low temperatures in Fig. [3J 

Conclusions: We have presented a theory for half 
metal-superconductor heterostructures and have inves- 
tigated the Josephson coupling through a half-metallic 
layer with a thickness of several coherence lengths. The 
Josephson coupling is induced by triplet pairing correla- 
tions in the superconductor. These triplet correlations 
are coupled to the singlet superconducting order param- 
eter via a spin-rotation effect, which occurs when quasi- 
particles in the superconductor are reflected from a spin- 
polarized medium. We have performed self-consistent nu- 
merical calculations for this problem, and found a low- 
temperature anomaly in the temperature behavior of the 
critical Josephson current. This anomaly is a robust fea- 
ture, which is not very sensitive to parameter variations. 
We discuss the Andreev excitation spectrum in the half 
metallic region, and explain the temperature variation of 
the Josephson current in terms of this spectrum. 
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